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ABSTRACT 


Lax-Wendroff  methods  for  hyperbolic  systems  have  two  characteristics 
which  are  sometimes  troublesome.  They  are  sonetimes  too  dissipative — 
they  may  smooth  the  solution  excessively- -and  their  dissipative  behavior 
does  not  affect  all  modes  of  the  solution  equally.  Both  of  these 
difficulties  can  be  remedied  by  adding  properly  chosen  accretive  terms. 

We  develop  modifications  of  the  Lax-Wendroff  method  which  equilibrate 
the  dissipativity  over  the  fundamental  modes  of  the  solution  and  allow 
the  magnitude  of  the  dissipation  to  be  controlled.  We  show  that  these 
methods  are  stable  for  the  mixed  initial  boundary  value  problem  and 


develop  analogous  formulations  for  the  two-step  Lax-Wendroff  and 
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MacCormack  methods. 
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1. 


Introduction 


We  shall  consider  approximations  for  hyperbolic  systems  of  partial 
differential  equations.  We  begin  by  considering  the  Cauchy  problem  for 

(l.l)  u,  = Au  , -®  < X < » , t > 0 

Xi  X 

with  initial  data 


vl.2)  u(x,0)  = U^(x)  , -00  < X < oo  , 

where  u is  a vector  of  length  n and  A is  an  nxn  matrix 
with  real  eigenvalues.  We  assume  that  A has  a complete  set  of 
eigenvectors  and  can  therefore  be  transformed  to  diagonal  form.  We 
vill  denote  the  eigenvalues  of  A by  If  A is  a function 

of  x,t  we  assume  that  this  transformation  can  be  done  smoothly. 

We  shall  discuss  the  well-known  Lax-Wendroff  method,  several  of 
its  variants,  and  modifications  thereof.  Discussions  of  these  methods 
and  modifications  of  them  which  improve  their  phase  errors  and  stability 
regions  have  been  carried  out  by  Turkel  [9],  Gottlieb  and  Turkel  [2], 
and  Eilon,  Gottlieb  and  Zwas  [1].  We  are  going  to  discuss  the 
dissipative  properties  of  these  approximations  and  modifications  of 
them  which  improve  their  dissipative  properties.  We  will  also  comment 
on  the  combination  of  our  modifications  with  some  of  those  of  the 
previously  mentioned  authors.  Recently,  Turkel  flo]  has  discussed 


I 


i 


a hybrid  leap-frog — Lax  Wendroff  method  which  is  less 
dissipative  than  the  Lax -Wendroff  method  but,  like  the  Lax -Wendroff 
method,  has  dissipation  which  varies  radically  for  the  various  modes 
of  the  solution. 

To  make  this  more  precise  we  now  introduce  some  notation  and 
definitions  which  are  discussed  in  Richtmyer  and  Morton  [ 8 ] and 
Kreiss  and  Oliger  [6], 

In  order  to  approximate  (l.l)  we  introduce  a grid  function 
v^(t)  = v(x^,t),  x^  = Vh,  h > 0,  V = 0,  + 1,  + 2,...  and 
t = 0,  k,  2k,..., k > 0.  We  write  our  approximations  in  the  form 


(1.3) 


v^(t+k) 


I Q V (t-jk) 
j=0 


where 


E and  Ev^ 

f =-00 


V 


V+1  ' 


Associated  with  (I.3)  is  the  characteristic  equation 


(1.4) 


Cr(p  = (h^'*’^ 


E Q.H'^)cp  = 0 
j=o 


where  Q,  = E A e | = (uh. 
by  letting  v(x,t)  = h ' e^^(p(cu). 
called  the  symbol  of  (I.3). 


Equation  (l.4)  is  obtained  from  (l.3) 
G,  as  defined  by  (I.4),  is  often 


Definition  1.1.  The  approximation  (I.3)  is  said  to  be  accurate  of  order 
^*^1^^^  for  solutions  u of  (l.l)  if  there  is  a function  C(t)  which 


3 


r 


k ■ 
) 


r-- 


is  bounded  on  every  finite  interval  [0,T]  such  that,  for  all 
sufficiently  small  h and  k, 


(1.5)  j|u(x,t+k)  - ^ Q.u(x,t-jk)(|  < kC(t)  (h  ^+k^) 


j=0 


If  <1]_  = Qg  ~ ‘I  approximation  (I.3)  is  accurate  of  order 

q- 

We  can  now  introduce  PCreiss'  definition  of  dissipativity  [6,  8]. 


Definition  1.2.  If  the  solutions  of  (I.5)  satisfy 


(1.6) 


|2r 


I K . I <1  - ?=)  I 5 1 for  0 < III  < 7T 
J 


for  some  5 > 0 and  natural  number  r,  then  the  approximation  (I.3) 
is  said  to  be  dissipative  of  order  2r. 

One  important  consequence  of  dissipativity  is  the  theorem  of  Kreiss 
and  Parlett  [6,  8]. 


Theorem  1.1.  Let  (I.3)  be  accurate  of  order  2m-2  or  2m-l  and 
dissipative  of  order  2m,  then  (I.3)  is  strictly  stable. 


Recall  that  strict  stability  implies  that  the  l^-nom  of  the  approx- 


imation, as  a function  of  t,  does  not  grow  faster  than  the  I^-norm 
of  the  solution. 


The  Lax-Wendroff  approximation  for  (l.l)  is 


(1.7) 


v^(t^-k)  = v^(t)  + kAD^v^(t)  + Y A^D+D_v^(t) 


1+ 


where 


D^VvU)  = (2h)-^(v^^^(t)-v^_^(t)) 


D^Vy(t)  = h'  (v^+;j_(t)-v^(t)) 


This  method  is  well  known  to  be  accurate  of  order  2 and  dissipative 

of  order  1 if  >,  = k/h  satisfies  0 < ?s  max  ln.j  < 1.  The 

0 

characteristic  equation  (I.5)  can  be  written  as 


(1.8)  K,  = 1 + iAn-  sin  | - 2A^n^  sin^(,?/2)  , j = 1,2,. ..,n 
JO  J 


(1.9) 


\k.\^  =1  - 4A^n^(l-7v^a^)  sin^l/2  si  - m.a^ 
J J J J 


where  we  let  a = sin  i/2  and  m.  = ^ . (l-l^^u  . ) • This  equality 

J J J 

yields  an  inequality  of  the  form  (1.6)  if  0 < 1 max  l^.l  < 1 but  it  is 

J ^ 

more  convenient  to  leave  it  in  this  form  for  our  purposes.  Note  that  the 

m.  are  functions  of  the  n-  and  that  the  "amount  of  dissipation",  the 
J <3 

amount  Ih.1  differs  from  1,  is  dependent  on  the  eigenvalues  of  A.  The 
J 

dissipation  is  greatest  for  intermediate  values  of  |u^l  and  least 


for  the  smallest  and  greatest  values  of 


See  Fig.  1 where  we 


plot  |h.|  as  a function  of  p.  = >4.  for  several  values  of  The 

J J J 

fact  that  the  dissipation  vanishes  for  all  a if  =0  is  often 

troublesome  in  nonlinear  calculations.  The  onset  of  nonlinear 
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w 


I 


» 

\ 
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instability  is  often  attributed  to  these  facts,  this  has  been 

discussed  by  Richtmyer  and  Morton  [ 8 ] beginning  on  page  334,  and  by 

Turkel  [ 9 ] • It  is  also  often  true  that  the  fastest  moving  waves 

associated  with  the  largest  values  of  ||ij|  are  often  of  minimal  interest, 

contain  the  greatest  observational  error,  and  are  contaminated  by  the 

most  computational  error.  In  such  situations  it  would  be  appropriate 

to  dissipate  these  modes  most  rapidly  rather  than  least  rapidly. 

Further,  there  is  generally  no  reason  why  intermediate  modes  should  be 

singled  out  for  most  rapid  destruction.  It  is  this  aspect  of  Lax- 

wendroff  methods  that  we  shall  discuss.  Our  subsequent  modifications 

will  be  directed  to  the  equilibration  of  dissipativity  over  the  modes 

lvalues  of  . ) or,  alternatively,  to  produce  decay  of  |k.|  as  a 
J J 

function  of  lu^.  1. 

Before  proceeding  we  shall  introduce  the  two-step  Lax-Wendroff 
method  and  the  MacCormack  method  which  are  a'lso  dissipative  and  behave 
similarly.  The  two-step  Lax-Wendroff  method  can  be  written  as 


+ V, 


V-l/2 


(t) 


(l.lO) 


Vylt+k)  = v^(t)  + kAD^l|W^ 
where  l/2  Vl/2  ^ ^ ' 

Our  introduction  of  indices  of  the  form  V + l/2  deviates  from  our 
earlier  definition  of  the  grid  function  but  the  meaning  should  be  clear. 
The  characteristic  equation  for  (l.lO^  is  just  (1.8)  in  the  linear 
case  we  are  considering  and  the  equality  (1.9)  holds  in  this  case  too. 

In  this  situation  (l.lO)  is  simply  a rearrangement  of  (1.7^.  However, 


i 


6 


oui-  modifications  of  (l.io)  will  not  be  the  same  as  those  of  (1.7)  since 
we  may  want  to  take  advantage  of  the  separate  steps  and  implement  them 
in  a two-step  manner.  The  MacCormack  method  can  be  written 


(1.11) 


% = Vv(t)  + kAD^v^(t) 


V (t)  + V,,  , 

v^(t+k)  = 2 + I AD_v^ 


where  the  are  intermediate  values.  The  characteristic  equation 


can  be  written 


(1.12) 


K = 1/2  + T + 2iAp.V3^T  + 2Au.a^T 

<J  J J 


where  T = 1/2  + a - Aq.a^  . 

J 0 

The  H can  again  be  seen  to  satisfy  (I.8)  and  (1.9).  Again,  in  this 
simple  situation,  this  is  a rearrangement  of  (1.7)  but  our  modifications 
will  again  be  different  since  we  will  implement  them  in  a two-step 
manner. 


7, 


2.  The  Modified  Methods  and  Their  Properties 

A nondissipative  method  like  the  leap-frog  method, 

Vy(t+k)  = v^(t-k)  + 2kAD^v^(t)  , 

can  be  modified  to  yield  a dissipative  method  by  adding  dissipative 
terms,  e.g., 

P 

v^(t+k)  = v^(t-k)  + 2kAD^v^(t)  " ^ v^(t-k) 

is  accurate  of  order  2 and  dissipative  of  order  for 

0<€<1,  |k|<l-e[6].  The  eigenvalues  of  the  symbol  of  this 

method  satisfy 

Uj  1^  = 1 - € sin^(|-) 

so  the  amount  of  dissipation,  the  magnitude  of  6 in  (1.6),  can  be 
controlled  by  varying  e. 

We  can  similarly  modify  the  Lax-Wendroff  methods  by  adding  terms-- 
we  cam  reduce  the  amount  of  dissipation  by  adding  accretive  terms. 

Such  terms  must  be  of  the  order  of  the  truncation  error  of  the  method 
so  that  they  do  not  constitute  a modification  of  the  differential 
equation  and  do  not  affect  the  order  of  accuracy  and  the  rate  of 
convergence  of  the  method. 

We  first  consider  a modification  of  the  Lax-Wendroff  method  (1.7). 
Let  and  be  arbitrary  matrices  which  are  diaigonalizable 

by  the  same  transformation  which  diagonalizes  A.  We  consider 


J 

:i 

( 

i 


8 


i- 


iL 


where  the  m *s  are  the  diagonal  entries  which  result  when  M is 
J -t 


transformed  to  diagonal  form  with  A. 

We  now  define  the  phase  error  per  time  step,  E,  of  an  approxi- 
mation as  (see  [ 9 ]): 


E = (approximations  phase  speed  - solutions  phase  speed)  x k 


Then  it  follows  that  [ 9 ] 


(2.1)-) 


E = arctan  ( (im  G)(Re  g) 


The  components  of  E = diag  (e^,...,e^)  for  our  diagonalized  system  are 


(2.5)  e^.  = ^ 0(§' 


We  now  consider  two  specific  modifications.  Let  = 0. 

If  we  take 


(2.6) 


2M^  = [-€l  + 4A^A^(I-A^A^)] 


t-hen  (2.3)  becomes 


(2.7) 


where 


2 2 

= -‘  +■  l4-p.(l-p.)  and  p.  = Au . 

u «3  0 0 


lo 
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We  have  tiius  cancelled  out  the  . dependence  of  the  coefficient 
of  cc^.  We  demonstrate  the  effect  of  this  in  Figures  2a-2d  where 
we  plot  |k.1  as  a function  of  (3.  for  several  values  of  | and  e. 
For  smaller  values  of  e the  dissipation  is  reduced  and  nearly  constant 
for  a considerably  larger  neighborhood  of  p = 0.  For  larger  values  of 
e the  dissipation  does  not  vanish  for  all  i in  the  neighborhood  of 


P = 0- 


We  next  consider  a modification  which  introduces  a quadratic  decay 


in  !k.|  as  a function  of  u..  We  take 
J J 


(2.8) 


2M^  = [-5>v^A^  + ] . 


Equation  \2.3)  is  now 


12.9) 


I |2  , 2 2 2 1+1  1+ 

U.l  = 1 - ]u 


where 


We  demonstrate  the  effect  of  this  modification  in  Figures  3a-3d 

where  we  plot  |k.|  as  a function  of  p.  for  several  values  of  5 
0 J 

and  & as  before. 

The  stability  of  the  modified  methods  given  in  (2.6)  and  1.2.9) 
follows  from  our  general  results  in  Gection  3-  We  will  state  the 
particular  form  those  results  take  for  these  methods  here.  The 
modification  (2.6)  is  strictly  stable  if 


f .. 

L 


(2.10) 


0 < ^ < i+(i-ip  r) 

m ax 


where  p is  the  value  of  p.  = A^.  with  largest  magnitude.  The 
max  j j 


modification  (2.9)  is  strictly  stable  if 


(2.11) 


0 < 6 < 


i^(i-lp  r) 

^ Tno  V ' 


We  could  have  introduced  a term  -8AA  in  (2.9)  instead  of  the  term 
-&A^A^  and  this  wovild  result  in  a linear  decay  of  |h.|^  with  respect 

t) 

4 

to  p.  in  the  a term.  However,  the  stability  analysis  corresponding 
J 

to  (2.11)  would  yield  0 < 6p.  < 4(l-p^)  for  j = l,2,...,n  so 

0 J 

that  6 would  have  to  be  chosen  to  agree  in  sign  with  g.,  i.e., 

J 

51  would  need  to  be  replaced  by  a matrix  which,  after  transformation 
to  diagonal  form,  would  yield  a matrix  with  entries  +5  in  appro- 
priate places  to  match  up  with  the  ji..  If  all  the  |j.  are  of  the 

J J 

same  sign  this  poses  no  problem;  however,  if  they  are  of  both  signs 

this  does  not  yield  a practical  procedure  in  general. 

We  now  note  that  the  first  term  of  our  expression  (2.5)  for  the 

phase  error  is  not  affected  by  our  choice  of  M^,  so  and 

can  be  chosen  as  in  Gottlieb  and  Turkel  [ 2 ] to  reduce  the  phase  error 

to  0(5^).  However,  these  modifications  do  affect  our  expression  for 

given  in  (2.j)  and  the  results  (2.I0)  and  (2.II)  no  longer  hold. 

J 

(2.3)  must  be  reexamined  when  nonzero  and  are  used  to 

establish  the  stability  bovinds  for  €,5,  |p  |,  etc. 

md,x 

We  now  consider  modifications  of  the  two-step  Lax-Wendroff  method. 
We  begin  by  noting  that  the  modifications  of  the  Lax-Wendroff  method 


which  we  have  already  discussed  can  all  be  used  in  the  second  step  of 
the  two-step  method  as  given  in  (l.io).  Since  the  symbol  is  unchanged 


12 


all  o\ax  previous  results  hold  for  such  modifications.  We  next  consider 
modifications  to  the  first  step  of  (l.lo)  as  given  by 


^ 2 ^ ^ 2 ^2^^^^  - IT  V2\(^) 


(2.12) 


v^(t+k)  = v^(t)  + kADj|)v^  + ^ M2hV(|)v^(t) 


The  symbol  for  this  method  is 


= 1 + 2XAia[-/l-a^ ' + i(AAa  + ^)]  + , 


(2.13) 


/V  2 2 2 k U 

Re(G)  = 1 - 2A  A^a  - M^Ma^  + f^a  , 


Im(G)  = 2MaVl-a^  , 


(2.14) 


1 . 


The  phase  error  is  now  given  by 


(2.15) 


We  again  note  that  e,.  is  not  changed  through  terms  of  order  5' 


p 

f 


Hni 


and.  that  the  phase  error  modification  schemes  of  Gottlieb  and  Turkel 
C8U1  again  be  applied  with  our  modifications  in  a straightforward 
manner. 

If  we  choose 


(2.16) 


2M^  = + 6AA 


and 


Mg  = 


= 0 we  have 


(2.17) 

with 


Notice  that  in  this  case  we  have  a quadratic  decay  of  |h.|  with  p. 
in  the  a term.  Thus,  this  method  is  similar  to  our  earlier  modifi- 


cation (2.8).  In  fact,  (2.17)  is  identical  to  12.9).  Our  remarks 
following  (2.9)  about  linear  decay  modifications  also  apply  here. 

The  stability  of  the  modification  (2.I6)  again  follows  from  our 
general  results  in  Section  3.  This  method  is  stable  if  the  inequality 
(2.11)  is  satisfied. 

We  now  t\irn  to  modifications  of  the  MacCormack  scheme.  For  the 
same  reason  as  before  for  the  two-step  Lax-Wendroff  method,  we  only 
consider  modifications  to  the  first  step  of  (l.ll).  Modifications  to 
the  second  step  will  be  the  same  as  those  for  Lax-Wendroff.  We  consider 
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r 


(2.18) 


= [I  + kAD^  + M^h^D^D_]v^i,t)  , 


Vy(t)  + V 

v^(t+k)  = — — + 2 ^ 


The  symbol  for  this  method  is 


(2.19) 


G = I + T + 2i>Aa  \Jl-a  T 2AAa^T 


where 


T = i + iXAVl-a^  - 'hha  - UiM^o^  \/l-a'  + , 


Re(G)  = 1 + a^(-2>pA^)+a^ll4-M^  + 8AAM^)  , 


Im 


(g)  = a'J  1-a  [2AA  + a^(-UM^)]  , 


and 


(2.20) 


IkJ^  = 1 + ^ bm  ]a^ 

J 0 J 

" [ ♦ I6m^j Vj  ^ I6\j 

+ + 64m^jXuj  ]a® 


The  phase  error  is  now  given  by 


(2.21) 


L '' 

■> 


Notice  that  the  o(|^)  term  of  the  phase  error  is  affected  by  our 
choice  of  this  time.  Simultaneous  modifications  to  improve  the 

phase  error  are  not  so  easy  to  carry  out  for  this  method. 
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w 


If  we  choose 

(2.22)  8M^  = - el 

I 

we  have  j 

U . = 1 - ea^  + O(a^)  I 

^ ! 

j 

5 

We  illustrate  the  effect  of  this  modification  in  Figures  Ua-hd  where  i 

we  again  plot  |h.|  as  a function  of  g.  for  selected  values  of  | and  e.  i 

We  can  also  choose  ' 

i 

(2.23)  8M^  = - &A^A^ 

to  introduce  a quadratic  decay  so  that 

= 1 - + o(a^) 

J J 

I 

We  have  also  plotted  this  in  Figures  5a-5d. 

The  situation  here  is  unlike  that  arising  from  our  previous  modifi- 
cations. Previously  the  stability  limit  for  increased  as  f 

decreased  to  zero  and,  in  fact,  approached  1 as  e tended  to  zero.  ! 

Here  the  stability  limit  tends  to  zero  as  e tends  to  zero.  From  the 
plots  in  Figures  l+a-Ud  it  seems  that  we  should  use  e somewhere  between 
0.1  and  2 in  order  to  have  a reasonable  stability  limit  for  |p  |.  The 

TT130C 

same  comment  applies  to  the  MacCormack  scheme  with  quadratic  decay. 

The  case  shown  in  Fig.  lib  seems  to  achieve  our  goal  of  equilibration  : 

quite  well  but  increases  the  dissipation. 

In  summary  we  remark  that  modifications  of  the  second  step  analogous 
to  those  used  with  Lax-Wendroff  are  much  more  successful  for  both  the  two- 
step  Lax-Wendroff  and  the  MacCormack  methods. 
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First  of  all,  note  that  all  the  unmodified  methods  ve  have  considered, 
namely,  the  one-step  Lux-Wendroff , the  tvo-step  Lax-Wendroff  and  the 
MacCormack  schemes,  have  the  same  ^ . and  are  all  dissipative  of  order 


t (and  therefore  strictly  stable  by  Theorem  l.l)  if  I ^ 

2 

each  of  our  modified  schemes,  -we  can  write  the  modified  k.  in  the 

J 

following  form: 


kJ  = 1 - a"^FO.,in^.,m2.,m3.,a2) 


By  Theorem  1.1,  our  modified  schemes  will  be  strictly  stable  if,  for  each  j. 


'Fi^yXa^yXn^y'm.^yC^)  > 0 for  a^e[0,l] 


Since  max  (5*2)  will  imply  (1.6). 

0 < I?  I < 7T  I sin 

Usually,  the  are  functions  of  Hence  (3.2)  defines  a 

stability  limit  for  B..  It  is  difficult  to  determine  this  stability 

u 

limit  for  general  functions  m^j(Bj)*  However,  for  specific  functions 
m^j(Bj),  this  stability  limit  can  always  be  determined  n\merically  as  in 
our  plots  of  |h. !• 

Consider  the  modified  one-step  Lax-Wendroff  scheme  (2.l)  with 
Mo  = M3  = 0.  The  following  theorem  gives  the  conditions  which  the 
will  have  to  satisfy  in  order  to  guarantee  stability. 


Theorem  3.1.  Consider  the  scheme  (2.1)  with  Mp  = M3  = 0. 


then  the  resulting  scheme  is  dissipative  of  order  1+  and  hence  strictly 
stable  by  Theorem  1.1. 

Proof;  Condition  (3.2)  implies  that  the  following  has  to  be  satisfied: 


F = -m^  .a^  + m (Up^a^-2)  +•  Up^(l-p^)  > 0 


Since  F is  a quadratic  in  m^^  with  negative  leading  coefficients, 
(3.^)  will  be  satisfied  if  (see  Fig.  6),  for  every  j, 

(i)  F has  real  roots,  say  mj^(o;^)  and  mg(o;^)  for  the  larger 
and  the  smaller  root,  respectively  (considered  as  functions  of  a),  and 


(3.3)  (ii) 


nic(a^)  HI  (a^)  for  all  a^e[0,l] 


= [-1  + -t-'Vl  + \!a  , 


(3-6) 


" 2p5'1-p5’ 


mc,(.a^)  = [ -1  2p^a‘"  -Tl  + - 1+cr^p^  Ma  , 

V J J 


= undefined 


if 

2 

a 

^ 0 

if 

2 

a 

= 0 

if 

2 

a 

^ 0 

if 

2 

a 

= 0 

Hence  the  roots  will  be  real  if 


Up  ? p 
1 Ua  p . - Ua  p . > 0 
J J 


or  if 

or  if 
13.7) 


for  all  j and  for  all  a e[0,l] 


for  all  j and  for  all  a €[0,l] 


2 1 

P < min  — 5 — = 1 for  all  j , 

a < 1 Ua  il-a  ) 


i.e.,  if  |p  |<1- 
' ' max  ' 


Note  that  (3-7)  is  just  the  unmodified  stability  limit.  Now,  by 
straightforward  differentiation,  we  can  show  that 


dm^  {cF ) dm  icF ) „ 

i,3.8)  — — >0  and  — — >0  for  a. -[0,1]  if  Ip  |<1 

dcF  ~ d«2  - 


Hence,  (3*5)  will  be  satisfied  if 


-2  il-p  . ) = m^.  U ) < m^  . < m (.O) 

Jo  IJ  b 


2 2 ' 
2p.a-p., 


for  all  j 


which  is  (3.3  ).  C] 


1 


Corollary  5.1:  The  modification  [2.6)  is  strictly  stable  if  1,2.10^  is 

— “ 

satisfied. 

Proof:  By  the  previous  theorem,  the  following  has  to  be  satisfied: 


-2(l-p^)  < - I + 2p^(l-p^)  < 2p^(l-p^l  for  all  j 

J ^ U J J J 


which  reduces  to 

(.5.9)  4(l-p*!^  > e > 0 for  all  j 

J 

Note  that  (3*9)  automatically  implies  |p  t < 1.  Also,  \2.10)  easily 

ind>x 

follows  frcm  (3.9).  Q 

Corollary  3-2:  The  modification  12.9)  is  strictly  stable  if  (,2.1l)  is 

satisfied. 

Proof:  Follows  immediately  from  Theorem  3«1‘  D 

Next,  consider  the  two-step  Lax-Wendroff  modification  12.12)  with 

= 0.  The  F we  obtain  in  this  case  is  very  similar  to  that  which 

we  obtained  for  the  Lax-Wendroff  method,  (3.^+).  We  only  need  to  replace 

ra,  . in  (3.*+)  by  m p.  to  obtain  the  correct  F.  Thus,  we  easily 
IJ  bj  3 

obtain  the  following  theorem. 


Tiieorem  3.2:  Consider  the  scheme  [2.12)  with  = 0.  If 

Ip. 


max 


< 1 


(3.10) 


and 


-2p^(l-3^)  < m .p  < 2(l-p^)  for  all 
J J J 


J ^ 


then  the  resulting  scheme  is  dissipative  of  order  k and  hence  strictly 
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w 


stable  by  Theorem  1.1. 

Corollary  3»3:  The  modification  (2.16)  is  strictly  stable  if  12.11)  is 

satisfied. 

A stability  analysis  is  more  difficult  for  modifications  of  the 

first  step  of  the  tiacCormack  method.  No  condition  analogous  to  the 

condition  (3.8)  holds  in  this  case.  We  have  not  been  able  to  obtain 

clean  conditions  like  those  in  Theorems  3*1  and  3 >2  for  the  MacCormack 

scheme.  However,  given  a specific  one  can  easily  determine  the 

stability  interval  by  examining  |h.1  as  a function  of  p.. 

J J 

We  finally  remark  that  the  previous  sufficient  conditions  for 
stability  are  all  necessary  if  we  allow  equality.  This  easily  follows 
from  the  von  Neumann  necessary  condition  [6,  8]. 
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U.  The  Initial  Bo\andary  Value  Problem 

We  now  consider  the  problem  of  approximating  equation  (l.l)  on  a 
bounded  x-  interval  a < x < b.  Since  the  modified  methods  we  have 
discussed  have  larger  stencils  (involve  more  neighboring  points)  than 
their  unmodified  counterparts  there  are  more  points  at  the  ends  of  the 
interval  [a,b]  where  these  approximations  cannot  be  used  than  there 
are  with  the  original  methods.  However,  this  problem  is  easily  avoided. 

Stable  approximations  for  the  initial  boundary-value  problem  for 
the  Lax-Wendroff , two-step  Lax-Wendroff , and  MacCormack  methods  are 
discussed  by  Gustaffson,  et.  al.  [4  ] and  by  Gottlieb  and  Turkel  [5  ]. 

We  will  base  our  methods  on  these. 


Assumptions 

We  assume  that  boundary  conditions  are  given  at  the  points  a 
and  b which  yield  a well  posed  problem  for  (l.l),  see  Kreiss  [ 5 ]• 

We  further  assume  that  stable  approximations  for  this  problem  are  known 
for  the  underlying  method  that  we  are  modifying,  see  [5,  4]  for 
candidates , and  finally  that  the  mesh  ratio  A and  modification 
parameters  e,  6,  etc.,  are  chosen  so  that  both  the  modified  and  un- 
modified methods  are  stable  for  the  Cauchy  problem. 

We  now  form  our  approximations  for  the  initial  boundary -value 
problem  by  coupling  the  \mmodified  method  with  its  stable  boundary 
conditions  to  the  modified  method  in  the  neighborhood  of  the  boundary 
points  a and  b in  the  manner  discussed  by  Oliger  [ 7 ] • 

The  Methods 

We  will  use  the  desired  modified  method  at  all  interior  net  points 


2? 


where  it  can  be  used,  we  then  d p the  modification  iset 
and  use  tne  underlying  method  at  all  those  points  in  the  neighborhood 
of  a and  b vdcere  it  can  be  used  iat  most  one  or  two  points  at  each 
end),  the  remaining  points  lonly  a and  b)  are  then  treated  using 
the  stable  boundary  approximation. 


Theorem  4.1.  The  methods  proposed  above  are  stable  in  the  sense  of 
Gustaffson,  et.  al.  [ 4 ] (definition  3-3)  if  our  assumptions  hold. 


Proof;  Tliis  result  follows  immediately  from  Tlieorem  2.b  of  [Y  ].  □ 
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